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Experimental Measurement of Mixed State Geometric Phase by Quantum 

Interferometry using NMR 

Arindam Ghosh and Anil KumaiQ 

NMR Quantum Computation and Quantum Information Group 
Department of Physics and NMR Research Centre, 
Indian Institute of Science, Bangalore-560012, India 

Geometric phase has been proposed as one of the promising methodologies to perform 
fault tolerant quantum computations. However, since decoherence plays a crucial role in such 
studies, understanding of mixed state geometric phase has become important. While mixed 
state geometric phase was first introduced mathematically by Uhlmann, recently Sjoqvist et 
al. [Phys. Rev. Lett. 85(14), 2845 (2000)] have described the mixed state geometric phase in 
the context of quantum interference and shown theoretically that the visibility as well as the 
shift of the interference pattern are functions of geometric phase and the purity of the mixed 
state. Here we report the first experimental study of the dependence of interference visibility and 
shift of the interference pattern on the mixed state geometric phase by Nuclear Magnetic Resonance. 

Key Words : Mixed state geometric phase, Interferometry, NMR. 



I. INTRODUCTION 



When a quantum system undergoes a unitary evolution and comes back to its initial state it 

hi , 

acquires a phase. The acquired phase can be of two types; the dynamic phase which depends on 
the evolution Hamiltonian and the geometric phase, which depends only on the evolution path 
of the quantum system in the projective Hilbert space For a two level quantum system 

(spin-|), the projective Hilbert space is a sphere and the geometric phase depends on the geodesical 
solid angle subtended at the center of the sphere by the path of evolution of the state vector. The 
concept of geometric phase first came in the adiabatic context P but later Aharonov et al. 
gave a non-adiabatic generalization of the theory of geometric phase . In the adiabatic approach 
the state vector is parallel transported adiabatically to ensure that the system always remains in 
one of the eigenstates (assuming that the system initially is prepared in one of the eigenstates) of 
the instantaneous Hamiltonian during the evolution. In the non-adiabatic approach the system 
is changed abruptly and the system comes back to its initial state through different intermediate 
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states. From the total phase acquired by the quantum system, the dynamical phase is eliminated 
by various methods in order to experimentally measure the geometric phase. In magnetic resonance 



events this i8 achieved by a spin echo Q. The theorv of phase of a pnre q uan t u m 

system or pure state geometric phase is well understood and has been demonstrated experimentally 



systems such as Nuclear Magnetic Resonance (NMR) single and two 



by various experimenta 
photon interferometry ?| 



Recently it has been proposed that fault tolerant quantum computation can be per 
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using geometric phase as it depends on the path and not on the speed of the evolution 
perform computation using geometric phase, it is necessary to understand the relation between 
geometric phase and decoherence. As decoherence or relaxation leads a system from a pure state to 
a mixed state, an understanding of the mixed state geometric phase is needed. In 1986 Uhlmann 
mathematically introduced the concept of mixed state geometric phase [10| . In this paper Uhlmann 
has taken a large system in pure state and a part or a subsystem in mixed state and pointed out 
the unitary evolution in which the subsystem is transported in a maximally parallel manner [lo| . 
Recently Sjoqvist et al. have provided a new description for mixed state geometric phase in terms 
of quantum interferometry 11]. In a quantum interferometer a quantum system undergoes a series 
of unitary evolutions, after which the probability of finding the system in one of its eigenstates 
becomes an oscillatory function of some control parameter. The oscillation pattern of the probability 
resembles the well known optical interference pattern. According to Sjoqvist et al. the shift of 
interference pattern is a function of the geometric phase acquired by the quantum system during 
the nnitary evo.nrions, as well as the purity of the initial interna! state (snch as, polarization of a 
photon) of the quantum systems involved in the interferometric operation [11| . The geometric phase 
can be directly measured from the shift of the interference pattern. Mixed state geometric phase 
has been experimentally demonstrated by Du et al. using NMR and by Ericsson et al. 
using single photon interferometry. Du et al. have experimentally demonstrated the mixed state 
geometric phase by measuring the relative phase change of an auxiliary spin. In the present work 
we measure the shift of the interference pattern in the Sjoqvist's interferometry model and show 
that the shift is same as the theoretically predicted geometric phase as a function of mixed state 
purity. We also demonstrate the effect of mixed state geometric phase on the interference visibility. 
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II. THEORY 
Quantum Interference 

Let us consider the Sjoqvist's interferometry model as shown in Fig^ Photons entering the 
interferometer along a horizontal path are split into two perpendicular paths by a beam splitter 
(BSi). On the horizontal path the photons are globally phase shifted, whereas on the other path 
the internal states of the photons, say the polarization, undergo a unitary evolution U. The photons 
are reflected by two mirrors (Mi,M 2 ) and the two paths meet again at another beam splitter (-BS2). 
A detector detects the photons coming only along the horizontal path. The detected intensity shows 
an interference pattern as a function of the phase shift. 

As a photon can take one of the two possible paths and in each path it can have one of the two 
possible polarizations, so the Hilbert space of the combined "path-internal state" system becomes 
2 2 x 2 2 . In NMR, the above interferometry model can be simulated using two coupled spin-| nuclei, 
which have the Hilbert space of identical dimension. One qubit represents the path qubit and the 
other qubit, termed as spin qubit, represents the internal state. 



The equivalent quantum circuit of the Sjoqvist's interferometry model is shown in Fig|2] |x3 ] . 
The eigenstates |0) and |1) of the path qubit represent the two paths, horizontal and vertical respec- 
tively. The path qubit is prepared in the pure state |0)(0| at the beginning of the interferometry 
operation. The beam splitter is represented by a Hadamard gate given by, 



1 1 
1 -1 



As the phase shifter and the unitary operator 'U' are path specific, they are represented by two 
controlled operations, together given by 
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® u + 
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(2) 



Depending upon the state of the first qubit, the operator Uq either applies U on the second qubit 
or phase shifts the first qubit. The mirrors in FigQ represent NOT gates, given by, 



U 



M 



1 

1 



(3) 



The input state of the combined "path-spin" system can be written as [14], 

p in = (\0)(0\) P ®Pl (4) 

where P stands for path and p$ is the density matrix corresponding to the initial state of the spin 
qubit which can be either pure or mixed. The initial density matrix pi n is transformed into the final 
density matrix p out as, 

Pout = U H U M U C U H p m U^U C ] U M ^U H ] . (5) 



Substituting the matrix forms of the operators given by Eq.(^|SJ) we obtain 11 | 



Pout 



1 1 
1 1 



1 -1 
-1 1 



1 1 

-1 -1 



1 -1 
1 -1 



<8> Upo 



(6) 



The detected signal in the horizontal path (|0) eigenstate of path qubit) is given by the trace of 
the reduced density matrix of the spin qubit corresponding to the |0) state of the path qubit. The 
output intensity is given by jl^ . 

= l - (1 + \Tr s (Up s ) \cos [ X - arg Tr s (Up s )]) 

= \ + u cos lx-4>}), 



(7) 



where the amplitude of oscillation v = \Tr$ (Upy) \ is called the visibility of interference and the 
shift (f> = arg Trs (Upo) depends on the unitary operator 'U 7 acting on the spin qubit density 
matrix p$. 



A mixed state can be thought of a mixture of several pure states incoherently weighted by 
their respective probabilities. Therefore, the interference pattern given by Eq|7| takes the following 
form for a mixed input spin statejll|. 



k \ k 



WkVk cos [x - 4> k ] 



where the index k denotes individual pure states with probabilities w^. The above equation can be 
written in the form of Eq0 as 



I oc 1 + v cos(x ~ 0) 
by defining mixed state phase shift <fi and visibility v as 



arg 



E 



w k u k e 



(9) 



(10) 
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Geometric phase and parallel transport condition 

The parallel transport condition for any state vector \if>(t)) is given by, 

<^)|#)>=0, 



(12) 



which means that the phase does not change when \tp(t)) evolves to \ip(t + 5t)) for infinitesimal St. 
When a mixed state given by the density matrix p m (t) evolves under a unitary operator Aft), the 
condition given in EqH21 leads to 



Tr 



Pm 

(t)A(t)A\t) 



0. 



(13) 



This condition although necessary is not sufficient to determine the unitary operator A(t) as it can 
determine A(t) only up to N phase factors, where N is the dimension of the Hilbert space. The N 
phase factors can be determined from the conditions jl5]. 



{k{t) A{t) A\t) k{t)) = 



k = 1,2,3. 



N, 



(14) 



where \k{t)) are the orthonormal eigenstates of p m {t). The unitary operator A(t), obtained by solving 
the above conditions, parallel transports the mixed state density matrix p m so that the dynamical 
phase becomes identically zero. The geometric phase 7 „ acquired by a mixed state when the state 
evolves under A(t) along a curve T, is given by 



7<?[r] = arg Tr[p m A(t)} = arg 



w k v k e 



(15) 



where e l/3fc is the geometric phase associated with the k th pure state. The expression for the 
geometric phase given by Eqll5lis similar to the expression for the interferometric phase shift given 
by EqEl and therefore the interferometric phase shift directly gives the geometric phase of the spin 
qubit. 



In the present work we consider the mixed state of a spin-| particle. The density operator 
of a spin-| particle can be in general written as, 

1 



Pr, 



[1 + f.a) 



(16) 



where the length 'r' of the Bloch vector r, is equal to one for pure states, less than one for mixed 
states and remains unchanged during unitary evolution of the state. The components of a are the 
Pauli matrices, a = [<j x , a y , a z ]. p m represents a mixture of two of its eigenvectors with eigenvalues 



|(l±r) 



Let us consider that the Bloch vector r for a mixed state (r<l) traces out a cyclic curve in the Bloch 
sphere which subtends a geodesically closed solid angle of Q. During the process the two eigenstates 
of the density operator with eigenvalues \ (1 ± r) acquire geometric phase =F^ respectively |l6[]. The 
quantity ^ k w k v k e %< ^ k fEqlTUlandfTT]). using the fact that v k = \ for cyclic evolution becomes, 

^liW** = i(l-r)^ + ^(l + r) e ^ 

k 

= cos(—) — irsin(—). (17) 

Using the expression given by EqEfrhe shift of interference pattern fEq JTU|) and the visibility fEq lTT]) 
for mixed state can be respectively written as, 

arctan | r tan ( — 

and 



cos 2 (— ) + r 2 sin 2 (— ) . (19) 

In the present paper, we have experimentally measured the above frequency shift and the visibility 
using a two qubit system, by NMR. The frequency shift directly gives the geometric phase of the 
spin qubit. 



III. EXPERIMENTAL PROCEDURE 



Experiments were performed on Carbon-13 enriched CHCl?, dissolved in CDCI3. The C and 
X H nuclei form a two qubit system with a J-coupling of 209 Hz. 1 H and 13 C respectively are used as 
the path and spin qubits. The spin-lattice (Ti) relaxation times of 13 C and 1 H at room temperature 
were measured as 21s and 16s respectively, and the spin-spin (T 2 ) relaxation times were measured to 
be 0.29s and 3.4s respectively. All the experiments were performed using a Bruker DRX 500 MHz 
(11.2 Tesla) NMR spectrometer where the resonance frequencies for 13 C and l H are 125.76 MHz and 
500.13 MHz respectively. The pulse programme is given in FigOl The pulse programme contains 
several parts which are described below: 
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Creation of Pseudo-pure state (PPS) 

The "path- spin" system is first prepared in a pseudo-pure state using the method of spatial aver- 



aging 



17| . The pulse sequence is as follows, 

where the superscript H or C identifies the spin (proton or carbon respectively) on which the r.f. 
pulse is applied and the subscript xor y determines the phase of the pulse. Jqh is the J-coupling and 
G z indicates a ^-gradient which destroys all coherences (x and y magnetizations) and retains only 
longitudinal magnetization (z magnetization component). At the end of this sequence the system is 
prepared in the 1 00) pseudo pure statei 

Creation of Mixed state 

After preparing the 1 00) PPS, an a degree pulse is applied on the carbon spin followed by a z- 
gradient. In the Bloch representation it creates a mixed state vector whose length r = cos a [r < 1 for 
0° < a < 90°], where the value of r determines the purity of the state. The above pulse programme 
can be written as, 

(af x -Gz. (21) 

The Interferometer 

I — I 

The Hadamard gate (Beam Splitter) is implemented by the sequence — (tt)^ [18] . The im- 

portant operation of the interferometry part is the controlled operation Uq- Uc contains a controlled 
phase shift gate applied on the path qubit and a controlled U operation acting on the spin qubit. 

Controlled Phase shift: 

The controlled phase shift gate in the present context is different from the conventional two-qubit 
gate. Here both the controlling and the target qubits are the path qubit. The path qubit is phase 
shifted by \ when it is in the |0) (horizontal path) state. The output of the phase shift gate U x can 
be written as, 

U x \00) = e*|00); U x \0l) = e *|01); U X \W) = |10); U x \ll) = |11). 
The pulse sequence for U x is |li 



GO - (*>? - (?) (22) 
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Controlled U: 

In the present case U is a geometric phase shift operator. The controlled geometric phase shift 
operator is implemented by evolving the spin qubit in a cyclic path (the 'Slice Circuit') in the Bloch 
sphere as shown in FigEJ when the path qubit is in state |1) (vertical path). This is achieved by 
pulsing only on the |10) — |11) subsystem. Two transition selective it pulses are applied on the 



1 10) — 1 11) transition of 13 C with phases differing by (it + <j)) Q|. They cause the Bloch vector to 
flip along one path, I\ and come back to its initial orientation along a different path, T 2 . The loop 
(ri,T 2 ) subtends a geodesically closed solid angle of Q = 20 P|. A spin echo sequence has been 
incorporated to eliminate the dynamical phase. The pulse sequence is given by, 

/ / s|10)-|ll) / x|10)-|ll) / /oq n 

where r is the total time of the two transition selective pulses. The second (it)? pulse in the above 
sequence restores the sign of the l H magnetization. 

The mirror (NOT gate) is implemented by a (it) pulse. It converts the state |0) (horizontal 
path) to state |1) (vertical path) and vice- versa. The sequence (■§) — {ft)x f° r t ne Hadamard gate 
is repeated after the mirror, in order to implement the second beam splitter (-65*2). 



Measurement 



At the end of the interferometric operations, both the qubits were detected after applying a z- 
gradient and a reading it/ 2 pulse on the detection qubit. The diagonal part of the density matrix 
was then tomographed using the line intensities normalized to the respective equilibrium spectra 



IV. RESULTS 



The intensity of the signal detected only in the horizontal path is proportional to the total pop- 
ulation of the 1 00) and 1 01) levels as these two energy levels correspond to the |0) state (horizontal 
path) of the path qubit. For each \ the final density matrix (diagonal part only) was tomographed 
and the sum of the 1 00) and 1 01) populations was plotted as the intensity. Data were collected at 
37 equidistant values of \ ranging from -360° to 360° to obtain the full interference pattern. Fig|S] 
shows the l H and 13 C spectra for different x values for pure initial state of the spin qubit. The 
three low intensity lines in the carbon spectra arise from the natural abundant carbon coupled to 
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deuterium in the solvent CDCI3. The intensities (sum of 1 00) and |01) populations) calculated from 
the normalized spectral line intensities are plotted as a function of phase shift \. Fig0a) shows the 
pattern when U (the geometric phase shift operator) was not applied. As expected no shift in the 
pattern from x — was observed. Whereas FigJS^b) shows the pattern corresponding to ft — 180°. 
A shift of -90° was observed as expected according to EqUHlfor pure state (r = 1). In each plot the 
solid line represents the expected theoretical curve. 

The 'shift - geometric phase' relationship (Ea ll8ll 

For mixed input state of the spin qubit, the pattern shifts from x — for non-zero geometric phase. 
The amount of shift is a function of both the purity of the input state as well as the geometric phase 
of the spin qubit. Fig shows the dependence of interferometric shift on the geometric phase and 
the purity of mixed state for ft = 60° (0a), ft = 90° (0b) and ft = 120° (0c). For a particular value 
of ft and r, experiments have been performed for ten equidistant points of x hi the range [-90°, 0°]. 
The data were fitted with a function ^(v, <p) — v cos(x — 4>), to calculate the shift. The shift is zero 
for r = and the shift is — ^ for r = 1 in all the three cases. The solid line in each plot represents 
the theoretical curve. Spectra corresponding to a = 0°, 30°, 50°, 70° and 90° have been shown for x 
= 30° (0a), x = 40° (0b) and x = 60° (0c). 

The 'visibility - geometric phase' relationship (Ea ll9B 

The visibility of interference or the amplitude of oscillation is given by the difference between the 
maximum and the minimum intensities in the interference pattern. The visibility was measured for 
different purity 'r' of the spin qubit state. Fig0 shows the visibility as a function of V for ft = 120° 
(0a), ft = 180° (0b) and ft = 360° (0c). ft = 360° makes the visibility independent of r, while for ft 
= 180° the visibility changes linearly with r. The experimental data matches the expected behavior 
(solid line) given by EalTTIl Spectra corresponding to a = 0°,30°,60° and 90° (for ft = 180°, a = 89° 
was applied instead of 90° as the shift according to EqlTHl becomes undefined for ft = 180° and r = 
0) have been shown adjacent to each plot. While recording the spectra the value of x was chosen 
according to the shift of pattern given by Eqll81 All the experiments for ft = 360° was performed 
at x = 0° an d for ft = 180°, at x = -90°. For ft = 120°, x was chosen same as the shift predicted 
by EqUHl 
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V. CONCLUSION 

The study of mixed state geometric phase has become important ever since geometric phase was 
proposed as a possible method of performing fault tolerant quantum computing. The pure state 
geometric phase is well understood and well studied by various experimental methods. Here we 
have reported the first experimental measurement of mixed state geometric phase directly from the 
shift of a quantum interference pattern. We have experimentally measured the visibility and the 
shift of the interference pattern as a function of the purity of the input mixed state which agree 
with the theoretically expected results. This study shows that NMR interferometry is one of the 
possible experimental methods to measure geometric phase of a pure as well as a mixed state. Future 
directions include studies of non-cyclic geometric phase j^] and applications of geometric phase in 
fault tolerant quantum computations. 
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FIGURE CAPTIONS 

Figure 1. The Sjoqvist's interferometry model ^j]. BSi and BS2 are two beam splitters which split 
a photon beam in two perpendicular paths. Mi and M 2 are two mirrors. The unitary operation 
l U' acts on the internal state of the photons traveling along the vertical path while the photons 
traveling along the horizontal path are phase shifted by x by the phase shifter. The detector detects 
the photons coming along the horizontal path only. 



'_igure 2. The quantum equivalent circuit of Sjoqvist's interferometry model using two qubits 



14j . One qubit corresponds the path while the other corresponds the spin (internal state). The 



single qubit gates used are Hadamard (Uh, EqJTJ) and NOT (Um, Eq|3J). The control operation Uc 
consists of a control phase shift gate acting on the first qubit and a control 'U' operator acting on 
the second qubit (Eq|2J) • 



Figure 3. The pulse programme to perform NMR interferometry. The black and empty boxes 
represent tt/2 and ir pulses respectively while the gray boxes represent pulses with flip angle given 
on the top. The phase of a pulse is given at top of the pulse. The gray Gaussian shaped pulses are 
transition selective soft it pulses applied on the 1 10} — |11) Carbon transition. Gz is the z-Gradient 
pulses applied to kill all transverse magnetizations. J is the coupling between l H and 13 C and is 
equal to 209 Hz in the present case. The pulse programme can be divided in four major parts, 
namely the PPS, preparation of mixed state, interferometer and the measurement, which are 
discussed in detail in the text. The last ir/2 pulses (shown by broken lines) on each spin is the 
detection pulse and is applied at the detected spin one at a time. 



Figure 4. The slice circuit for cyclic evolution of a state vector [for pure state \r\ = 1 and 
for mixed states \r\ < 1] in the Bloch sphere jlflj . In the present case, the 1 10} — 1 1 1 } subsystem is 
transported through a closed loop using two transition selective tt pulses applied at phases (9) and 
(# + 7r + 0) on the 1 10) — 1 11} Carbon transition. The first tt pulse takes the second qubit from state |0) 
to state |1) along the path Ti, while the second tt pulse brings the state back from |1) to |0) alonga 
different path T 2 . The loop [Fi,r 2 ] subtends a solid angle f2 = 20 at the center of the Bloch sphere 5(. 



Figure 5. The l H and 13 C spectra obtained for various values of \i starting from pseudo 
pure initial state of the spin qubit. The additional low intensity lines in the carbon spectra arise 



13 

from the natural abundant carbon coupled to deuterium in the solvent CDCI3. The sum of 1 00) and 
1 01) populations, calculated from the normalized line intensities have been plotted as intensities 
against 37 equidistant values of the phase shift x m the range [-360°, 360°]. (a) Q = 0°. No shift 
in pattern was observed, (b) Q = 180°. A shift of -90° was observed. The solid line in each plot 
represents the theoretically expected value [EqJ7]. 

Figure 6. The shift in interference pattern as a function of purity of mixed state for Q = 
60° (a), Q = 90° (b) and Q = 120° (c). In each plot the shift is plotted for 6 different values of r 
(cosa) between and 1. 13 C and l H spectra are shown for 5 different pairs of a and \. The solid 
line in each plot represents the theoretically expected curve according to EqlTKl 

Figure 7. The dependence of interference visibility on mixed state purity for Q = 120° (a), 
Q = 180° (b) and Q = 360° (c). In each plot the visibility is plotted for 9 different values of r (cosa) 
ranging from to 1 (for Q = 180°, a was varied from 0° to 89° instead of 90° in order to avoid the 
singularity of Eq ll8|) . 13 C and l H spectra are shown for 4 different values of a. The solid line in 
each plot represents the theoretically expected curve according to EqfT9l 
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